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ABSTRACT 



We present results from an Ha monitoring campaign of the Be X-ray binary systems 
HDE 245770 = A 0535+26 and X Per. We use the Ha equivalent widths together with adopted 
values of the Be star effective temperature, disk inclination, and disk outer boundary to deter- 
mine the half-maximum emission radius of the disk as a function of time. The observations of 
HDE 245770 document the rapid spectral variability that apparently accompanied the regenera- 
tion of a new circumstellar disk. This disk grew rapidly during the years 1998 - 2000, but then 
slowed in growth in subsequent years. The outer disk radius is probably truncated by resonances 
between the disk gas and neutron star orbital periods. Two recent X-ray outbursts appear to 
coincide with the largest disk half-maximum emission radius attained over the last decade. Our 
observations of X Per indicate that its circumstellar disk has recently grown to near record pro- 
portions, and concurrently the system has dramatically increased in X-ray flux, presumably the 
result of enhanced mass accretion from the disk. We find that the Ha half-maximum emission 
radius of the disk surrounding X Per reached a size about six times larger than the stellar radius, 
a value, however, that is well below the minimum separation between the Be star and neutron 
star. We suggest that spiral arms excited by tidal interaction at periastron may help lift disk gas 
out to radii where accretion by the neutron star companion becomes more effective. 
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1. Introduction 

Be X-ray binaries (BeXRBs) consist of rapidly 
rotat i ng B -stars with neutron star companions 
llCoel l2000h . The B-star primaries lose mass 
into an outflowing circumstellar disk, and if the 
disk reaches a radius comparable to the perias- 
tron separation, then disk gas accreted by the 
neutron star can power a significant (and often 
trans i ent) X-ray so urce (jOkazaki fc Negueruela 



20011: lOkazaki et al.l l2002). The disk gas is de- 
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tected through observations of emission lines, 
an optical/infrared excess, and a net linear po- 
lariza tion caused by starlig ht scattered in the 
disk (jPorter fe Riviniud 120031 ). Be star disks 
are inherently time variable and can develop 
and disappear on tim e scales o f years to decades 



Hubert fc Floquet 



(lUnderhill fc Doazanl Il982t 
19981 ). Thus, we expect that the X-ray ac- 
cretion fluxes will vary on similar timescales. 
Here we present a joint study of the disk Ha 
emission flux and X-ray flux variations in two 
of the best studied BeXRBs, HDE 245770 = 
A 0535+26 and X Persei. We discuss similar co- 
variations of Ha emission and X-ray flux in com- 
panion papers on t he black hole bi nary Cyg X- 
1 = HDE 226868 (|Gies et all 120031). the super- 
giant binary LS I +65 010 (jGrundstrom et al 



2007a), and the BeXRB microquasar LS I +61 303 



(jGrundstrom et alj|2007bh 



The remarkable Be star HDE 245770 (and X- 
ray counterpart A 0535+26) is one of the pro- 
totype objects of the BeXRB class. The primary 
star of HDE 245770 has a classificatio n of 09.7 Hie 
Giovann elli fc Sabau Graziatil 119921 ) or B0 Hie 
Steele et al.l 19981) . but the luminosity class is 



still somewhat controversial since many of the lu- 
minosity criteria are based o n spectral lines tha t 
are affected by disk emission dWang fc Gieslll998l) . 
The orbital period fro m peaks in the X-ray flux is 
appro ximately 110 d (jHutchingsl Il984t ICoe et aT 
2006), and a comparison of the times of maxima 



between the first and most recent epochs suggests 
a period of P = 109.96 ± 0.05 d. The only reli- 
able estimate of the other orbital elements comes 
from X-ray timing of the puls ar spin period dur- 



ing three outbursts in 1993 (jFinger et alj 11994 ; 



Bildsten et al.l 119971 ) , and this solution indicates 
that the orbital eccentricity is e = 0.47 ± 0.02. 
The radial velocity curve of the Be star has an 
amplitude that is smaller than or comparable to 
the velocity fluctuations introduced by emission 
contamination and other factors, so only upper 
limits are available for the Be star's semiampli- 
tude and the system mass ratio, M 2 /Mi < 0.12 
(|Wang fc Gied fl998). The star has a long and rich 
history of X-ray, optical, and infrared observations 
dMotch et al.l 1991 ; Giovannelli fc Sabau Graziatil 



1992; Clark et al. 1998; Ncgueruela et al. 



Haigh. Coe. fc Fabregatl 12004 ICoe et al 



record of the light curve variations back to the year 
1898, and they show that after a quiescent stage 
in the first part of the century the star brightened 
by some 40% in the early 1970s. The circumstellar 
disk that appeared then provided the gas to power 
the X-ray source at the time coincident with the 
dawn of X-ray astronomy. 

The emission from the disk of HDE 245770 
in the form of optical and infrared continuum 
light and Balmer emission lines displays large 
temporal variations in the record since the 1970s 
( Clark et alll998tlLyuty fc Zaitseval2000HHaigh et al 



20041 ). The dominant timescales of variability are 
~ 1500 d and 103 d, where the latter is the beat 



perio d of the long and orbi tal periods (jHaigh et al 



20041 ). lHaigh et al.l (|2004l ) argue that the disk IR 



flux tends to hover around three different bright- 
ness levels, and they suggest that these corre- 
spond to distinct disk radii that are defined by 
resonances bet ween the disk gas and neutron star 
orbital periods. lOkazaki fc Negueruelal (|200ll ) de- 
veloped a model to predict the resonant disk trun- 
cation radii in BeXRBs, and lHaigh et al. I (|2004l ) 
found that the model predictio ns are consistent 

"~ ( 20061) 



Coe et al 



with the IR magnitude jumps 
extended this analysis with additional JHK pho- 
tometry and measurements of the Ha emission 
line . The Ha emission variations ( Clark et al 



1998; H aigh et al J 1 19991: iLvutv fc ZartsevaT koOO: 



Haigh et all 120041 : Icoe et all l2006h indicate that 
the disk almost totally disappeared in 1998 but 
rebounded to its former strength over the next 
few years. 

The binary system X Per (HD24534; B0 Ve, 



Lvubimkov et al.l 119971 ) is the brightest and per- 



haps most famous member of the BeXRB class. 
X Per has displayed time variable optical emis- 
sion lines since the earliest photographic spectro- 
grams w ere made at the be ginning of the last 
century (jCowlev et alJ [l972|). The long term 



Lvutv fc Zaitseval (|200oh present a remarkable 



cha nges in the s pectra l appe arance are r e viewe d 
by iRoche et al] (1 19931) and IClark et all (|200ll ). 
In a seminal paper, IClark et al.l (|200ll ) describe 
the photometric and spectroscopic variations ob- 
served over the period 1987 through 2001, and 
during this time the star made a remarkable trans- 
formation f r om B e to B and back again to Be. 
IClark et all (|200ll) discuss how these variations 
are related to the structural properties of the cir- 
cumstellar disk. The orbital elements were deter- 



2 



mined by iDelgado-Marti et al. (|200lh by careful 
timing observations of the X-ray pulsar (pulse 
period « 835 s) using the Rossi X-Ray Timing 
Explorer (RXTE) satellite. The orbital period 
is P = 250.3 d and the orbital eccentricity is 
e = 0.11. 

Both HDE 245770 and X Per are among 
some hundred objects regularly observed with 
the RXTE All-Sky Monitor (ASM) instrument 
( Levine et al. Il996l ). and here we address the is- 
sue of how their disk size variations are related to 
the observed X-ray fluxes. We present our recent 
observations of their Ha emission lines in §2. We 
then show in §3 how the Ha emission strength 
is related to the disk radius, and we use these 
relations in §4 to document how the disk radius 
variations relate to the y-band and X-ray light 
curves. The disk radii appear to attain limits 
that are probably related to reson ant truncation 
radii (jOkazaki fc Negueruelal 120011 ). and we dis- 
cuss in §5 how such limits affect the mass transfer 
process in these binaries. An appendix to the pa- 
per offers a reassessment of the orbital motion of 
X Per based upon radial velocity measurements of 
UV spectra from the International Ultraviolet Ex- 
plorer (IUE) satellite, from which we derive limits 
on the binary's mass ratio and inclination. 

2. Observations and Ha Variations 

We observed both HDE 245770 and X Per be- 
tween 1998 August and 2000 December with the 
Kitt Peak National Observatory 0.9 m coude feed 
telescope. We used two spectrograph arrange- 
ments to record the red spectral region around 
Ha with resolving powers of R = A/AA = 4100 
and 9500. The details about the spectra and their 
reduction for these runs are given in a compan- 
ion paper (jGrundstrom et al.ll2007af) . In addition 



to these primary runs, we obtained several more 
spectra during auxiliary runs in 2004 October and 
2006 October (using grating B in second order for 
a resolving power R — 9500 over the range 6470 - 
7140 A and 6433 - 7143 A, respectively). All the 
spectra record the Ha emission line and the He I 
A6678 feature. 

The general trends in the spectra are illustrated 
in Figures 1-3. We show in Figure 1 the set of 
spectra of HDE 245770 from the 1998 August - 
September run that coincided with the time of disk 



loss. The spectra are arranged with their continua 
set at the time of observation and scaled so that 
1 d along the y-axis corresponds to 10% of the 
continuum flux. We also show a model spectrum 
at the bottom of the figure that is derived from the 
grid of non-local thermodynamic equilibrium and 
line-blanketed model atmospheres calculated by 
Lanz fc Hubeml §003). This synthetic spectrum 



is based on the parameters of effective temper- 
ature T c ff = 28000 K, gravity logq = 3.3, and 
solar abundances (jGiovannelli fc Sabau Graziati 
[1992J). The spectrum was broadened by a sim- 
ple convolution with a rotational broadening 
function usin g a linear limb dark ening coeffi- 
cient of 0.24 (jWade fc Rucinski|[l985h and a pro- 
j ected rotational velocity of V sin i = 230 km s _1 



dGiovannelli fc Sabau Graziatill 992; Wa ng fc Gics 
1998HHaigh et al.ll2004 . We see that there is evi- 
dence of residual emission components in both the 
Ha and He I A6678 absorption lines at this time. 
The profiles also display significant night-to-night 
or faster variations. For example, there are well 
defined structures within the core of He I A6678 
that may result from photospheric nonradial pul- 
sations (see the ca se of the similar star £ Oph; 
Kambe et al.Ml997r ). and there are variations in 



the relative strengths of the violet and red emis- 
sion peaks in Ha. These suggest that mass loss 
processes leading to the development of the new 
disk were already underway at this time. Note 
that the model indicates that a line blend should 
be present for the C II A6578, 6582 doublet, which 
appears to be absent from the HDE 245770 spec- 
tra. These lines weaken at higher temperature, 
but an increase in temperature would result in a 
model He I A6678 profile weaker than observed. 
The absence of the C II featur es (and the weakness 
of the C IV profiles in the UV: IWang fc Giedll998h 
may point to a carbon deficiency caused by the 
presence of CNO-processed gas in the atmosphere 
of HDE 245770. 

The disk emission strength of HDE 245770 had 
increased dramatically by the time of the next runs 
in 1999. We show in Figure 2 the average spectrum 
from each of the runs with the continuum levels 
offset for clarity. The growth in the Ha emission 
is accompanied by the appearance of emission in 
He I A6678 that almost totally obscures the profile 
during the runs in 2000. Note how the separation 
between the violet and red peaks of Ha decreases 
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as the emission strengthens. 

We show in Figure 3 the average spectrum of 
X Per from each of the runs, and we see that Ha 
and He I A6678 were emission lines during this 
time. The overall morphology and strength of the 
lines agree well w ith the mor e exten ded time se- 
ries presented bv_ Clari^ira2j 2001 ) . Inspection 
of Figure 7 in IClark et al.l (|2001r ) shows that our 
observations were made while the disk was active 
and generally increasing in strength. 

We made a number of measurements to charac- 
terize the profile variations of Ha for both stars, 
and these are summarized in Table 1 (given in 
full in the electronic version). Column 1 lists 
the target name, column 2 gives the heliocen- 
tric Julian date of mid-exposure, and column 3 
shows the co rresponding orbital phase (from the 
ephemeris of iFineer et al.lll994l for H DE 245770 



and from iDelgado-Marti et al 



20011 for X Per). 



Column 4 lists the equivalent width determined 
by a numerical integration of the line flux over 
the range 6536 - 6590 A. We measured the 
radial velocity of the wings based upon a bi- 
se ctor position determined using the method 
of IShafter. Szkodv. fc Thorstensenl (|l986l ) . This 
method samples the line wings using oppositely 
signed Gaussian functions and determines the 
mid-point position between the wings by cross- 
correlating these Gaussians with the profile. We 
used Gaussian functions with FWHM = 200 
(100) km s" 1 at sample positions in the wings 
of ±320 (±190) km s" 1 for HDE 245770 (X Per), 
and these radial velocities are giv en in column 5 of 
Table 1. IZamanov et al. ( 1999h advocated mak- 
ing fits of double-peaked Ha profiles using Gaus- 
sian functions to match the violet V and red R 
peaks, and we have followed their approach here, 
although we caution that such a functional fit is 
poor in many instances. These double-Gaussian 
fits were restricted to the inner part of the profile 
(|AA| < 6 A for HDE 245770 and < 2 - 4 A for 
X Per) since the wings are much more extended 
that those of Gaussian functions. The remaining 
columns in Table 1 list the parameters for these 
fits: radial velocity of the V peak (col. 6), ra- 
dial velocity of the R peak (col. 7), ratio of the 
equivalent widths of the V and R components 
(col. 8), V/R peak intensity ratio (col. 9), FWHM 
for the V peak (col. 10), and FWHM for the R 
peak (col. 11). Note that Ha has only one peak in 



the 2006 observations of X Per, so we fit a single 
Gaussian to the profiles of these three spectra. 

We find no compelling evidence that these fit- 
ting parameters vary with orbital phase, although 
we caution that our orbital phase coverage is 
incomplete. In particular, the radial velocity 
measurements for the Ha wings and emission 
peaks do not appear follow the orbital motion ex- 
pected for the Be star (see lWang fc Gieslll998l for 
HDE 245770 and the Appendix for X Per). We 
suspect that these velocity measurements are dom- 
inated by fluctuations related to the azimuthal dis- 
tribution and non-Keplerian motion of the disk 
gas, and these variations confound attempts to 
measure the much smaller orbital motion in our 
data sets. The best strategy for measuring the or- 
bital radial velocities of these Be stars may be to 
wait until an episode of disk loss occurs and then 
embark on a program of high dispersion, high S/N 
spectroscopy of emission-free lines in the blue part 
of the spectrum. 

3. Disk Radius and Ha Strength 

The Ha emission forms primarily in the disk 
surrounding the Be star, and the total flux of the 
feature (measured as the line equivalent width) is 
closely rela ted to the size of the dis k. In a re- 
cent paper ( Grundstrom fc G~iesl 2006h we present 
calculations of the disk size and emission strength 
for an extensive range of stellar and disk proper- 
ties. The Be disk models are based upon a sim- 
ple parameterization of d isk p roperties outlined 
by iHummel fc Vranckenl ((2000), from which we 
derive a synthetic Ha line profile and the angu- 
lar distribution on the sky of the wavelength in- 
tegrated Ha flux. The model predictions about 
the emission strength and angular size of the disks 
agree with observational results for bright Be stars 
whose disks have been resolved through optical 
long baseline interferometry. 

The disk model assumes an axisymmetric gas 
structure with the gas density distribution set by 
a power law R~ m (where R is the distance from 
the rotational axis), an exponential decline above 
and below the disk, and a base density at the star's 
surface. The disk surface brightness is set by its 
temperature, and we assume an isothermal disk 
with a gas temperature of 0.6T o ff where T e g is the 
effective temperature of the Be star. The spatial 
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flux distribution is integrated along segments nor- 
mal to the major axis of the projected disk to ob- 
tain a collapsed emission sum as a function of po- 
sition along the major axis. We adopted a working 
definition of disk radius as the distance along the 
major axis from the star where the integrated Ha 
emission intensity of the disk declines to half of the 
peak value found at a position immediately adja- 
cent to the photosphere (the half-width at half- 
maximum intensity radius or HWHM radius for 
short). 

The functional relationship between Ha equiv- 
alent width and disk HWHM radius for a specific 
Be star depends upon the star's T e ff, the inclina- 
tion of the disk normal to the line of sight, and 
the adopted outer bounda ry for the disk radius 
(|Grundstrom fc Gied [2006h . We list in Table 2 
the values of these parameters that we adopt for 
HDE 245770 and X Per. The listed radii and 



masses are estimates from Okazaki fe Negueruela 



( 2001 ) and lLvubimkov et alj (ll997T) and the stellar 
effecti ve te mperatures are fromlGiovannelli fc Sabau 
(|l992h and lLvubimkov et all |l997j) for HDE 245770 
and X Per, respectively. We assume that the disk 
is co-planar with the orbit, and we adopt the mid- 
range estimates of inclination from Wang fc Giesl 
(|l998l ) and the Appendix for HDE 245770 and 
X Per, respectively. We set the outer boundary for 
the model disk equal to the Be star's Roche radius 
at apastron. The results are generally insensitive 
to this last assumption, since the Ha flux is small 
from the outer, optically thin portions of the disk. 
The predicted relationships between equivalent 
width W\ [Ha] and the ratio of the disk HWHM 
radius to the stellar radius Rd/R s are shown in 
Figure 4 for a large range in disk base density. 
The model equivalent width is referenced to the 
photospheric flux, and if the disk contributes a 
continuum flux fraction of e = F^/F^ then the 
observed equivalent width must be prorated by a 
factor of (1 + e) in order to find the ratio of the 
disk to star radius Rd/R s from Figure 4. 

There are several features of these simple mod- 
els that need to be considered in the derivation 
of a disk HWHM radius from the Ha equiva- 
lent width. First, if the disk temperature dif- 
fers from our assumed value of 0.6T c fj, then the 
predicted Ha equivalent width will change by a 
factor of approximately Tdisk/(0-6T e gf) for a given 
disk radius. This is due to the fact the disk gas 



source function varies almost linearly with tem- 
perature for Ha in hot Be stars like these two 
targets. Second, the equatorial gas density is as- 
sumed to follow a power law R~ m with m = 3.0, 
the value often derived in studies of the IR flux 
excess of Be stars. The numerical relationship be- 
tween Wa[Ho:] and Rd/Rs is relatively insensitive 
to th e selected value of m ( Grundstrom fc G~ies1 
2006), although the base density for any position 
along the curve depends critically on m. We cau- 
tion, however, that the actual disk density prop- 
erties may vary significantly from a power law 
in these BeXRBs. For example, if the star ex- 
periences an episode of increased disk mass loss, 
a density enhancement may appear to propagate 
outwards over time and produce increased emis- 
sion strength when the enhancement reaches the 
radius of the optically thick/thin boundary. Us- 
ing the disk radius derived from the relation in 
Figure 4, we would correctly determine that the 
HWHM radius had increased but we would err if 



GrazifflS assumed that the emission strengthening re- 
sults from an increase in density at all radii. 

Finally, it is possible the disk may be truncated 
by tidal or resonant effects at a radius smaller than 
our assumed outer boundary. If the disk density 
drops to zero in the outer optically thin zone, then 
only a small part of the total emission is lost and 
we would only slightly overestimate the emission 
flux for a given HWHM radius. Thus, the revised 
curves would appear slightly above those shown 
in Figure 4 (see an example of the variation with 
outer boundary in Fig. 1 of I Grundstrom fc Gied 



2006) and our derived HWHM radii will be slightly 
smaller than the actual values. On the other hand, 
if the disk density vanishes at a radius where the 
gas is optically thick, then the revised curves will 
reach an asymptotic limit slightly below the cut- 
off radius. The resulting functions would fall be- 
low the standard curves at the high emission end 
in Figure 4, and we would then tend to overes- 
timate the HWHM radii using the standard rela- 
tions. The former case probably applies to the two 
targets we analyze here, since the HWHM radii 
we find are generally small compared to the sus- 
pected truncation radii (we discuss this further in 
§4). Thus, the HWHM radii measurements we de- 
rive may slightly underestimate the actual values. 
Our analysis is insufficient to detect the presence 
of a truncation radius, but if one exists, its radius 
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will be larger than the HWHM radius. 

4. Disk Growth and X-ray Accretion Flux 

The mass transfer rate and subsequent accretion- 
driven X-ray flux will clearly depend on the 
changes in the radius of the Be star's disk that 
we can track through the variations in the Ha 
equivalent widths. Here we consider the long term 
variations in the Be star's disk as observed in the 
Ha emission flux and in the disk continuum flux 
(as seen in the observed F-band light curve). We 
then compare the derived disk radius variations 
with those observed in the X- ray light curves fro m 
the RTXE/ ASM instrument (jLevine et al.lll996f ). 



4.1. HDE 245770 



Lvutv fe Zait scva (2000) note that the F-band 



continuum of HDE 245770 may be much brighter 
than the stellar continuum alone (due to contin- 
uum light from the disk), and consequently the ob- 
served equivalent width referred to the combined 
continuum flux may underestimate the absolute 
Ha emission flux. In order to compare the Ha 
emission flux to a constant stellar continuum, we 
need to rcscale the equivalent width by the factor 
1 + e where e gives the ratio of the monochromatic 
disk-to-stellar c ontinuum flux at wavelen gths near 
Ha. Following [Lvutv fc Zaitseval (|2000h . we can 
estimate the rescaled flux by comparing the V 
magnitude at the time of the Ha measurement 
with the magnitude of the star alone, V = 9.50, 

by 

(1 + e ) W x = lO" - 4 ^" 9 - 50 ) W\. (I) 

Lvutv fc Zaitseval (|2000h present V magnitudes 
for dates contemporaneous with the Ha measure- 
ments through 1998, and we have used these to 
make the small corrections from equation (1). 
However, there are no published estimates of V 
for the succeeding years, so instead we relied on 
surrogate measur ements of the J magnitude from 
Coe et ahl (I2006T) in order to make these correc- 



Haigh et al.1 (|2004 (see their Fig. 1) show 



tions. 

how the infrared magnitude variations track those 
observed in the F-band for HD E 245770, and we 
used a series of photometr ic V ( Lvutv fc Zaitseval 
2000) and J observations ( Coe et all 20061) over a 
common time span to find a relation between these 
for HDE 245770, 



which has an empirical scatter of ±0.0 3 mag. We 
used t he observed J magnitudes from ICoe et aD 
(200G) with equation (2) to find V and transform 
the equivalent widths from 1999 onwards to an 
absolute scale. We then used Figure 4 to find esti- 
mates of the disk HWHM emission radius for the 
times of the Ha observations. 

The observed light curve and derived disk 
HWHM radii are plotted in the middle and 
lower panels of Figure 5. The disk HWHM 
radii are derived from o ur W\ data together wit h 



like measurements from iGi ovanncll i et al.l (11 999). 



Lvutv fc Zaitseval (l200 0h. IPiccioni et all (l2000al). 
Haigh et all (|2004h and lcoe et all (|2006l ). We see 



(V - 9.50) = 0.32( J — 8.67), 



(2) 



that the disk almost disappeared in 1998 and then 
started to grow rapidly in radius over the next two 
years. The expansion rate slowed considerably by 
2001 at a radius of Rd/Rs ~ 4. However, subse- 
quent observations show that the disk continued 
to grow slowly and reached a maximum radius of 
R d /R s w 5 by 2005 . 

ICoe et ah ( 20061 ) also present a time evolution 
diagram of disk radius (see their Fig. 7) that they 
derive from the velocity separation Ay of the 
violet and red peaks of the Ha profile ( Huand 
19721 ). Their method leads to somewhat smaller 
radii than our HWHM radii (f» 83% as large us- 
ing AV measurements derived from Table 1) be- 
cause the radius derived from AV is based upon a 
velocity sampling weighted by the brighter parts 
of the disk while the HWHM radius corresponds 
to an outer boundary between the bright and faint 
parts of the disk. The two depictions of disk evolu- 
tion are in general agreement except for estimates 
from 2005 August where Coe et al. determine a 
disk radius about twice as large as indicated in 
our Figure 5. We think this discrepancy is due to 
the fact Coe et al. relied on an extrapolation of a 
linear relationship between WApa] and AV^ to es- 
timate radius for the 2005 August data (for which 
no direct Ay measurements were available). If 
one adopts a st andard power law rela tion between 
these variables ( Zamanov et al. 200 ll ) rather than 
a linear fit (which tends to underestimate AV for 
strong emission), then the extrapolation leads to 
radii similar to those shown in Figure 5 for 2005 
August. 

We also show in the lower panel of Figure 5 
dotted lines that indicate disk HWHM radii de- 
rived from mean levels of Ha strength in the 



G 



past that were highlighted in the review by 



Lvutv fc Zaitseval (|2000h (and indicated in their 
Fig. 3). The lower two levels correspond to the 
two groupings of Ha strength observed between 
1987 and 1998 while the upper level shows the 
much larger disk HWHM radius associated with 
data from 1975 through 1981 when the Ha emis- 
sion was exceptiona lly strong (W\ ~ —26 A 
de Loore et al .1 119841 : V ps 8.9, iLyutv fc Zaitseval 



2000). Note that the deceleration in disk growth 
observed in 2001 occurred at the time when the 
disk HWHM radius had reached the lower cluster 
level. 

Haigh et all 1)20041 ) and lCoe et all (|2006h argue 



that the tendency for the disk emission fluxes to 
cluster at specified levels is related to the pres- 
ence of resonances between the disk gas and neu- 
tron star orbital periods th at tend to truncate the 
disk at specific disk radii ( Okazaki fc Negueruelal 
20011 ). These truncation radii are given by 



Rn = ( GM S 
Rs V 4^ 2 



1/3 



1 (P 



Rs 



orbit 



2/3 



r(M s /M ) 1 /3 



nV3R s /R Q 
(3) 

where n is the integer number of gas rotation peri- 
ods per neutron star orbit and r is a constant equal 
to 97 if the resonance occurs with the orbital pe- 
riod or 92 if the resonance is with the shorter beat 



period (103 d; lLarionov. Lvutv. fe Zaitseval 120011 ) 
presu mably caused by re tr ograde precession of the 
disk jHaigh et al.1 12004 ). lOkazaki fc Negueruelal 
(2001) predict that the important truncation radii 
for HDE 245770 will be those associated with the 
n = 4 and 5 resonances, and using their adopted 
parameters (Table 2) the resonance radii will oc- 
cur at Rd/R s — 6.6 and 5.7 for n — 4 and 5, 
respectively (for r = 92). The n = 5 truncation 
radius is indicated as a dashed line in the lower 
panel of Figure 5, and we see that the disk HWHM 
radius hovered slightly below this limit between 
2001 and 2006. This suggests that the rapid disk 
growth phase ended when the outer parts of the 
disk reached the n = 4 or 5 resonance truncation 
radius. 

The top panel of Figure 5 shows the mean X- 
ray fluxea_J binned in time slots equal to one or- 
bital period, and error bars indicate the standard 
deviation of the mean within the time bin. The 



3 http: / /xte. mit.edu 



slow expansion of the disk apparently led to the 
large X-ray outb urst observed near JD 2,453,526 
(|Coe et al.l 120061 ) and a second smaller outburst 
near JD 2,453,616. These outbursts occurred 
when the disk had reached its largest HWHM ra- 
dius over the duration of the RXTE/ ASM mis- 
sion. The binary semimajor axis is approximately 
17.9R S (Table 2) and the periastron separation 
is 9.5i? s (shown as the dot-dashed line in lower 
panel of Fig. 5). The mean Roche radius of the 
Be star will be approximately 5.7i? s at periastron, 
about the size of disk HWHM radius at the time 
of the outbursts. This suggests that the disk had 
grown to a size sufficient to permit active mass 
transfer at periastron around the epoch when the 
outbursts were observed. Furthermore, the huge 
emission strengths observed during 1976 - 1981 in- 
dicate that the disk may have attained a HWHM 
radius of Rd/R s = 9 then (top dotted line in the 
lower panel of Fig. 5), exceeding the Roche limit 
at periastron. 

We suspect that the mass transfer process is 
probably aided b y the strong tidal forc es that exist 
near periastron. lOkazaki et al.l (|2002l ) present hy- 
drodynamical simulations for BeXRBs that show 
how the periastron tides excite a two-arm spiral 
structure in the Be star disk, and the arm closest 
to the neutron star can lift material far beyond the 
nominal disk boundary to provide a source of gas 
accretion at phases beyond periastron. We show 
in Figure 6 the X-ray light curves bin ned in orbital 
phase according to the ephemeris of Finger et al.l 
(|19941 ) for the three energy bands observed by 
RXTE/ ASM. The time sample is restricted to qui- 
escent (non-outburst) dates. We see that the high 
energy flux attains a maximum near orbital phase 
0.3, i.e., well past periastron at phase 0.0. A simi- 
lar phase de lay in the emi s sion s trength of H/3 was 
detected by iMotch et al.l (|199lh . These observa- 
tions suggest that the disk does become more ex- 
tended for a period following periastron, leading to 
a delayed peak in mass transfer and X-ray emis- 
sion. 

4.2. X Persei 

The Ha emission equivalent widths of X Per 
observed in 2000 and 2004 and especially in 
2006 are the largest measured over the last few 
decades (Roche et al.lll993l : IPiccioni et al. 2000bl : 
Clark et al.l 2001 ). and they indicate that the disk 
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has grown significantly in size. Once again we 
need to rescale the equivalent widths to a con- 
stant stellar continuum flux level by considering 
the (1 + e) factor. We estimated this flux excess 
using the V-band light curve and the relation 



(1 + e) W x = io-°- 4 ^- 6 - 78 ) W x 



(4) 



where the base level magnitude of the star alone, 
V = 6.78, is that observed during the last disk- 
free phase ( Lvubimkov et al.lll997l ). The V-band 
light curve comes from observations by the mem- 
bers of the American Association of Variable Star 
Observers (AAVSO) that we transformed to stan- 
dard Johnson V magnit ude using contempora- 
neous photometry fro m Zamanov fc Zamanova 
(|l995l) . lEngin fc Yucel dl997j), an d photoelectric 
meas urements from the AAVSO ([Percy fe Bakosl 
20011 ). The transformed AAVSO measurements 



were binned into 25 d means to increase the S/N 
ratio. 

Our results are plotted as function of time in 
Figure 7. The lower panel shows the time evo- 
lution of the disk HWHM radius derived from 
our Ha equivalent width measurements and other 
publi s hed values dClark et al.l l200ll : iLiu fc Hand 
120011 IZamanovet al.l l200lh . The middle panel 
shows the time binned V-band light curve from the 
AAVSO observations. We see that X Per reached 
a maximum of V ~ 6.2 in 2000 and then rose 
again to V « 6.1 where it has remained to the 
present. The circumstellar disk of X Per is now 
brighter than it has been for the past few decades 



(jRoche et al.1 11993I ). The top panel shows the 



binned X-ray fluxes observed with RXTE/ ASM. 
We see that the first brightening episode corre- 
sponded to the onset of a slow increase in X-ray 
flux that eventually soared to a record high over 
the ASM observation period shortly after the sec- 
ond visual brightening. The X-ray flux continues 
to remain high during the current optically bright 
state. 

The apparent increase in X-ray flux (and im- 
plied gas accretion rate) that accompanied the 
disk expansion confirms that the X-ray source is 
powered by gas from the Be star disk. However, 
the largest disk HWHM radius shown in Figure 7 
is still much smaller than the separation between 
the Be star and neutron star. For example, the 
mean Roche radius at periastron is 34i?g (Ta- 
ble 2), which is a factor of five larger than the 



maximum disk HWHM radius (Rd/R s = 6.4 in 
Fig. 7). The maximum predicted disk size due 
to the action of the tidally-driven eccentric in- 
stability occurs at the 3 : 1 resonance radius be- 
tween the disk gas and neutron star orbital p eriods 



tween tne qisk gas ana n eutron star orbital p eriods 
(lOkazaki fc Neeueruelal I200I : IClark et al.l boOlh . 
and this radius occurs at Rd/R s = 31 (Table 2). 
Once again the observed maximum HWHM radius 
is smaller than the predicted truncation limit. 

However, tidal forces at periastron may excite 
a two-armed spiral in the disk, and the spiral arm 
facing the companion can lift disk g as out to close 



to th e vicinity of the neutron star (jOkazaki et al 
2002). If so, mass transfer will occur predomi- 
nantly after periastron as the gas in the arm ex- 
tension is accreted by the neutron star. We show 
in Figure 8 the X-ray light curves for X Per from 
fluxes recorded in the recent active state (HJD > 
2,452,200) th at are binned according to t he orbital 
ephemeris of Delgado-Martf et all (|200ll ) . We see 
that there is an orbital modulation (particularly in 
the high energy band) that was not obvi ous in the 
earlier, low state data ( Wen et aT1l2006h and that 
the X-ray maximum occurs about one quarter of 
a period after periastron. This phase of maximum 
is consistent with post-periastron accretion from 
extended disk gas in a spiral arm. 

5. Discussion 

The Ha observations presented here and else- 
where document the remarkable and continuing 
changes that occur in the mass loss from the Be 
stars in these two BcXRB systems. The gas that 
enters their circumstellar disks becomes a reservoir 
of fuel for accretion by the neutron star compan- 
ion, and the X-ray activity that accompanies the 
mass transfer increases dramatically as the disk ra- 
dius increases in size. The disk growth is expected 
to be limited by g ravitational interaction s with 
the neutron star (Okazaki fc Negueruelal 2001 
Okazaki et al. I l2002h . In low eccentricity bina- 



ries like X Per, the limiting radius occurs at the 
3:1 resonance radius by the tidally-driven eccen- 
tric instability, while in high eccentricity systems 
like HDE 245770, the neutron star will approach 
closer to the Be star at periastron and higher inte- 
ger resonances occurring at smaller radii will act 
to truncate the disk. 

Our observations appear to confirm these ex- 



S 



pectations in the case of HDE 245770. The largest 
disk HWHM radius we find from the Ha equiva- 
lent width is comparable to both the n = 5 res- 
onance radius and the mean Roche radius at the 
time of periastron. The historical maxima of Ha 
strength may imply that the disk HWHM radius 
can occasionally grow to even larger dimensions. 
We caution, however, that it may be possible for 
Be disks of extreme density to attain very large Ha 
emission fluxes without radial growth. For exam- 
ple, pressure broadening in high density disks can 
lead to an increased flux in the wings of Ha and 
thus a larger emission equivalent width. Whether 
such high density and high pressure disks could 
remain confined by gravitational forces is an open 
question. 

The case of X Per, on the other hand, shows 
that the mass transfer and X-ray accretion flux 
can increase even when the disk HWHM radius at- 
tains a size well below the critical value. We found 
that the recent X-ray flux increase occurred when 
the disk HWHM radius grew to about Rd/R s = 5, 
much less than the resonant truncation radius of 
R3/ R s = 31. There are two probable explanations 
for this difference. First, tidal forces at perias- 
tron promote the develo pment of a two-arm spi- 
ral structure in the disk (Okaza ki et al. 1 120021) and 
gas concentrated in the arm closest to the com- 
panion can be pushed outwards and accreted by 
the neutron star. Hydrodynamical models sug- 
gest that the gas accretion rate from the spiral 
arm will peak sometime after periastron, and we 
find that indeed the X-ray flux maximum lags pe- 
riastron by « 25% of the orbital period in both 
the BcXRBs examine d here and in the micro - 
quasar LS I +61 303 (|Grundstrom et al.ll2007bl ). 



Secondly, the disk probably extends beyond the 
HWHM radius of the Ha emission that we use 



to define disk radius (jGrundstrom fe; Giesl 120061 ) . 
The disk density is described in our model in terms 
of a power law with radial d istance from the star 
([Hummel fc Vranckenl [2000h . and the lower gas 
density at larger radius may be entirely adequate 
to power the accretion-driven X-ray flux. 

The observational record for both of these 
BeXRBs documents the growth of the disk from 
inside out. A comparison of the F-band light 
curves and Ha radii in both Figures 5 and 7 shows 
that the time evolution of the Ha emis sion lags be- 
hind that of the disk continuum flux. Clark et al 



(|200lh also noted this delay between the line and 
continuum fluxes, and they argue that the time 
lag is due to differences in spatial origin. The con- 
tinuum flux excess is probably formed in the inner 
part of the disk while the Ha emission is more op- 
tically thick and forms over a wider range of disk 
radius. Thus, an outwards propagating density 
enhancement would peak first in continuum light 
and later in the Ha flu x. The fate of t his ou tflow 
is unclear. Models by [Oka 

zaki et al.l (|2002h sug- 
gest that some of the mass falls back onto the Be 
star, some is accreted by the neutron star, and the 
remainder escapes from the system (presumably 
into a circumbinary region centered on the orbital 
plane). A search for such circumstellar structures 
might yield evidence of the ejected gas (for exam- 
ple, as a compact H II region surrounding X Per; 
iRevnolds et~al~ll2005l ). 
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A. Mass Ratio Limit from IUE Radial Velocities of X Per 



Several investigators have attempted to measure the orbital motion of X Per, but their results have 
generally been inconclusive because the lines are broad, shallow, and often marred by emission and because 
the semiamplitude is probably small (jrlutchingsl 11971 iRevnolds et~ai1 Il992t ISticklandl [l992). The ideal 
spectral range to search for orbital motion is in the ult raviolet because the disk flux contribution is relatively 
small in the UV ( Telting et al.|[l998h . Stickland (1992) used the spectra available at that time in the archive 
of the IUE to measure radial velocities, and he found that the velocity excursions were too small to measure 
orbital motion. Here we repeat the an alysis of the now larger set of IUE spectra using the complementary 
orbital elements from the pulsar orbit (|Delgado-Marti et al.l l 200lh . 

We obtained 43 high dispersion, short wavelength prime camera spectra of X Per from the IUE database 
maintained at the Multimission Archive at Space Telescop^]. These spectra were transformed to a uniform 
log A, heliocentric wavelength grid, normalized to a pseudo-contin uum, and excised of the main interstellar 
lines, and radial velocities were measured by cross-correlation ([Penny. Gies. fc Bagnuolol 1999h . A mean 
formed from the average of all the spectra was used as a spectral template, and its absolu te velocity was 
deter mined to be +0.8 km s -1 by cross-correlation with a similar spectrum of HD 34078 (jGies fc Bolton 
1986). We used the apparent shifts between each cross-correlation function and the ensemble mean of the 
functions to estimate the relative velocity shift for each spectrum, and then the final absolute velocities were 
set by adding the template velocity above. The results are listed in Table 3 (given in full in the electronic 
version) that gives the hel iocentric Julian date of mid-exposure, the orbital phase from the ephemeris of 
Delgado-Marti et al. I (|200lh . and the measured radial velocity. 

We then fit the IUE velocities using the non-linear, least-squares, orbital elements program of 



Morbev fc Brosterhud (Il974h. This was a constrained solution in which all the elements were set from 
the pulsar results (|Delgado-Marti et al. 200lh (with the longitude of periastron changed by 180°) with the 
exception of the systemic velocity 7 and the semiamplitude K. The formal solution for these parameters 
is 7 = 1.0 ± 0.9 km s _1 and K = 2.3 ± 1.4 km s _1 , and the observations and calculated radial velocity 
curve are illustrated in Figure 9. Clearly the amplitude of motion is small enough that this result is of 
marginal significance, but the velocities do place a useful upper limit on the semiamplitude. The dotted line 
in Figure 9 shows the fit with a semiamplitude of K + 2cr = 5.0 km s _1 , and this appears to be a reasonable 
upper limit consistent with the spread in the IUE velocities. 



The pulsar orbital elements (jDelgado-Marti et alj 120011 ) give a neutron star orbital semiamplitude of 



K = 39.8 ± 0.4 km s , so the ratio of the upper limit of the Be star semiamplitude to that for the pulsar 
yields an upper limit on the mass ratio, M2/M1 < 0.13. We have plotted this constraint in the mass 
plane diagram in Figure 10 together with lines of constant orbital inclination derived from the pulsar mass 
function. A lower li mit on the inclina ti on of i > 2 3° re sults from assuming that the star rotates slower 
than the critical rate (jClark et al.ll2001r ). lHarmaned (jl988T ) presents a summary of mass and radii data from 
eclipsing binary stars, and the mass range for a B0 V star is probably between ss 11M Q and 17M . The 
shaded reg ion in Figure 10 shows this probable Be star ran ge together with the observed range in neutron 
star mass ([van Kerkwijk. van Paradijs. fc Zuiderwij kill 9951 ). This region in the mass plane corresponds to 
an inclination range of i = 28° to 35°. 



7 http:/ /archive. stsci.edu/iuc/ 
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Table 1 
Ha Measurements 



Star 


Date 


Orbital 


W x 


V r (W) 


V r (V) 


V r (R) 


W X (V) I 




FWHM(V) 


FWHM(fl) 


Name 


(HJD-2,400,000) 


Phase 


(A) 


(km s- 1 ) 


(kms- 1 ) 


(km s- 1 ) 


W X (R) 


V/R 


(A) 


(A) 


HDE 245770 


51055.969 


0.103 


1.76 


18.4 














HDE 245770 


51056.989 


0.112 


1.96 


42.4 














HDE 245770 


51057.960 


0.121 


1.41 


32.4 














HDE 245770 


51058.950 


0.130 


2.16 


30.5 














HDE 245770 


51061.954 


0.157 


1.96 


38.0 














HDE 245770 


51061.993 


0.158 


2.07 


31.7 














HDE 245770 


51063.962 


0.176 


1.62 


21.1 














HDE 245770 


51065.981 


0.194 


1.05 


16.5 














HDE 245770 


51066.902 


0.202 


0.52 


18.0 














HDE 245770 


51066.976 


0.203 


0.88 


20.3 














HDE 245770 


51419.972 


0.403 


-3.11 


22.5 


-147.7 


165.5 


0.80 


0.96 


4.41 


5.34 


HDE 245770 


51419.995 


0.403 


—3.20 


23.1 


— 144.9 


165.6 


0.80 


0.95 


4.68 


5.52 


HDE 245770 


51425.987 


0.458 


-4.31 


8.6 


-158.3 


161.6 


0.95 


0.97 


5.12 


5.26 


HDE 245770 


51428.950 


0.485 


-4.06 


12.1 


-151.8 


165.0 


0.96 


0.98 


4.97 


5.10 


HDE 245770 


51464.958 


0.811 


—5.32 


4.5 


-157.6 


156.7 


0.98 


1.06 


4.84 


5.24 


HDE 245770 


51464.980 


0.811 


-5.14 


4.7 


-155.8 


158.8 


1.00 


1.05 


4.89 


5.14 


HDE 245770 


51465.936 


0.820 


-5.17 


3.1 


-156.5 


154.2 


0.99 


1.05 


4.75 


5.04 


HDE 245770 


51465.957 


0.820 


-5.04 


2.1 


-156.6 


153.9 


0.99 


1.03 


4.81 


5.01 


HDE 245770 


51466.906 


0.829 


-4.80 


-3.6 


-157.5 


149.1 


1.02 


1.04 


4.72 


4.78 


HDE 245770 


51466.928 


0.829 


-4.79 


-5.9 


-158.5 


149.3 


1.06 


1.03 


4.81 


4.69 


HDE 245770 


51467.988 


0.839 


-4.41 


-3.6 


-153.5 


155.5 


1.06 


1.05 


4.87 


4.85 


HDE 245770 


51468.952 


0.847 


-4.54 


-4.5 


-152.6 


159.5 


1.22 


1.06 


5.28 


4.61 


HDE 245770 


51468.974 


0.847 


-4.73 


-1.5 


-152.9 


155.8 


1.12 


1.06 


5.03 


4.75 


HDE 245770 


51469.916 


0.856 


-4.87 


-8.1 


-155.9 


152.1 


1.12 


1.06 


5.05 


4.76 


HDE 245770 


51491.904 


0.055 


-5.01 


6.0 


-148.6 


154.3 


1.01 


0.97 


4.96 


4.78 


HDE 245770 


51492.833 


0.064 


-5.20 


10.0 


— 148.0 


153.6 


0.93 


1.00 


4.74 


5.09 


HDE 245770 


51493.827 


0.073 


-5.25 


12.3 


-150.5 


156.9 


0.92 


0.98 


4.82 


5.10 


HDE 245770 


51494.869 


0.082 


-5.57 


12.8 


-147.5 


163.4 


0.95 


1.01 


4.76 


5.08 


HDE 245770 


51494.890 


0.082 


-5.36 


10.9 


-148.5 


156.6 


0.96 


1.01 


4.74 


4.98 


HDE 245770 


51495.944 


0.092 


-5.28 


6.7 


-148.6 


157.8 


1.01 


0.99 


4.92 


4.84 


HDE 245770 


51496.919 


0.101 


-5.20 


12.6 


-147.7 


161.0 


0.94 


1.00 


4.69 


4.98 


HDE 245770 


51497.852 


0.109 


-4.66 


17.8 


-144.8 


161.8 


0.92 


0.96 


4.76 


4.98 


HDE 245770 


51817.905 


0.011 


-8.03 


12.4 


-126.1 


133.5 


0.92 


0.96 


4.77 


5.00 


HDE 245770 


51818.908 


0.020 


-8.36 


6.2 


-120.7 


136.8 


1.04 


0.98 


5.11 


4.78 


HDE 245770 


51819.899 


0.029 


-8.00 


8.7 


-120.8 


136.4 


1.04 


1.00 


4.96 


4.77 


HDE 245770 


51821.884 


0.047 


-7.66 


6.6 


-121.6 


132.3 


0.99 


0.98 


4.83 


4.77 


HDE 245770 


51822.888 


0.056 


-7.83 


4.7 


-123.3 


133.3 


0.99 


0.98 


4.89 


4.84 


HDE 245770 


51823.830 


0.065 


-7.86 


11.9 


-119.7 


134.9 


0.98 


1.00 


4.84 


4.90 


HDE 245770 


51823.934 


0.066 


-8.22 


5.4 


-122.1 


134.8 


1.01 


1.01 


4.93 


4.92 


HDE 245770 


51824.855 


0.074 


-8.19 


13.4 


-121.5 


134.8 


0.94 


0.96 


4.90 


5.01 


HDE 245770 


51824.962 


0.075 


-7.86 


8.7 


-121.9 


131.6 


0.94 


0.99 


4.76 


5.02 


HDE 245770 


51830.913 


0.129 


-8.81 


13.2 


-115.1 


140.6 


1.09 


1.03 


5.15 


4.88 


HDE 245770 


51890.840 


0.672 


-7.21 


6.6 


-119.4 


130.1 


1.03 


1.01 


4.84 


4.76 


HDE 245770 


51890.861 


0.672 


-7.03 


7.3 


-120.6 


129.6 


1.00 


0.99 


4.79 


4.72 
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Star 


Date 


Orbital 


W x 


V r (W) 


V r (V) 


V r (R) 


W X (V) I 




FWHM(V) 


FWHM( 


Name 


(HJD-2,400,000) 


Phase 


(A) 


(km s _1 ) 


(km s _1 ) 


(km s _1 ) 


W\(R) 


V/R 


(A) 


(A) 


HDE 245770 


51892.806 


0.690 


-6.97 


6.1 


-119.8 


127.1 


0.98 


0.99 


4.79 


4.83 


HDE 245770 


51892.827 


0.690 


-6.99 


5.4 


-117.8 


130.3 


1.01 


0.99 


4.88 


4.76 


HDE 245770 


51893.812 


0.699 


-7.60 


6.5 


-117.2 


131.6 


1.04 


0.99 


4.96 


4.71 


HDE 245770 


51893.833 


0.699 


-7.55 


7.2 


-114.8 


134.4 


1.10 


0.99 


5.03 


4.56 


HDE 245770 


51894.799 


0.708 


-7.61 


5.7 


-117.9 


132.9 


1.08 


0.99 


5.08 


4.67 


HDE 245770 


51894.823 


0.708 


-7.65 


5.4 


-118.5 


132.5 


1.05 


0.99 


4.98 


4.70 


HDE 245770 


51895.834 


0.717 


-7.63 


7.5 


-118.2 


131.4 


1.04 


1.00 


5.00 


4.80 


HDE 245770 


51895.913 


0.718 


-7.53 


7.4 


-117.5 


131.8 


1.04 


1.00 


4.95 


4.77 


HDE 245770 


51896.834 


0.727 


-7.37 


7.5 


-118.7 


130.3 


1.00 


0.99 


4.84 


4.79 


HDE 245770 


51896.855 


0.727 


-7.46 


5.7 


-117.9 


132.1 


1.05 


1.00 


5.00 


4.75 


HDE 245770 


51897.832 


0.736 


-7.47 


9.2 


-119.2 


128.6 


0.98 


1.00 


4.81 


4.94 


HDE 245770 


51897.853 


0.736 


-7.57 


8.9 


-118.9 


129.9 


1.00 


1.00 


4.85 


4.86 


HDE 245770 


51898.841 


0.745 


-7.77 


8.0 


-118.8 


132.2 


1.03 


1.00 


4.89 


4.77 


HDE 245770 


51898.862 


0.745 


-7.78 


8.2 


-119.5 


132.0 


1.01 


1.01 


4.88 


4.85 


HDE 245770 


51899.838 


0.754 


-7.90 


8.3 


-118.8 


132.8 


1.02 


0.99 


4.96 


4.81 


HDE 245770 


51899.859 


0.754 


-7.86 


8.5 


-118.8 


134.2 


1.03 


0.99 


5.00 


4.80 


HDE 245770 


51900.831 


0.763 


-7.95 


6.1 


-117.8 


136.2 


1.11 


1.00 


5.18 


4.69 


HDE 245770 


51900.852 


0.763 


-8.00 


6.3 


-120.2 


135.0 


1.07 


1.00 


5.14 


4.83 


HDE 245770 


51901.817 


0.772 


-7.87 


8.3 


-115.8 


136.5 


1.11 


1.00 


5.18 


4.65 


HDE 245770 


51901.838 


0.772 


-7.79 


9.4 


-117.2 


135.0 


1.06 


0.99 


5.09 


4.77 


HDE 245770 


54019.981 


0.975 


-10.08 


0.5 


-119.4 


103.1 


0.82 


0.82 


4.92 


4.94 


HDE 245770 


54020.953 


0.984 


-10.21 


2.8 


-116.3 


105.8 


0.80 


0.84 


4.77 


4.99 


HDE 245770 


54021.957 


0.993 


-10.16 


3.1 


-113.7 


105.9 


0.84 


0.83 


4.88 


4.81 


X Per 


51055.896 


0.811 


-9.38 


4.5 


-73.1 


93.0 


0.79 


0.92 


2.71 


3.16 


X Per 


51056.966 


0.815 


-9.09 


4.6 


-72.3 


94.5 


0.77 


0.91 


2.71 


3.22 


X Per 


51057.948 


0.819 


-9.64 


4.4 


-71.2 


93.5 


0.82 


0.92 


2.74 


3.09 


X Per 


51058.938 


0.823 


-9.50 


4.0 


-70.4 


95.1 


0.85 


0.92 


2.81 


3.06 


X Per 


51061.986 


0.835 


-9.51 


0.6 


-67.6 


92.5 


0.99 


0.93 


2.92 


2.74 


X Per 


51062.915 


0.839 


-9.49 


0.8 


-69.1 


91.1 


0.93 


0.93 


2.84 


2.85 


X Per 


51063.954 


0.843 


-9.09 


1.2 


-68.9 


92.5 


0.90 


0.90 


2.85 


2.85 


X Per 


51065.936 


0.851 


-9.25 


0.7 


-67.2 


92.0 


1.00 


0.90 


2.96 


2.66 


X Per 


51066.937 


0.855 


-8.99 


-1.1 


-68.1 


90.4 


0.96 


0.90 


2.86 


2.70 


X Per 


51491.878 


0.553 


-10.05 


-15.1 


-116.4 


62.0 


0.71 


0.77 


3.65 


3.97 


X Per 


51492.818 


0.557 


-10.15 


-11.8 


-113.5 


64.2 


0.70 


0.77 


3.58 


3.94 


X Per 


51493.813 


0.561 


-10.24 


-10.6 


-110.2 


67.6 


0.76 


0.80 


3.72 


3.91 


X Per 


51494.845 


0.565 


-10.31 


-12.5 


-111.4 


65.9 


0.77 


0.80 


3.74 


3.91 


X Per 


51494.853 


0.565 


-10.27 


-12.4 


-112.3 


65.5 


0.75 


0.80 


3.68 


3.92 


X Per 


51495.873 


0.569 


-10.30 


-11.7 


-109.0 


67.8 


0.78 


0.82 


3.67 


3.85 


X Per 


51496.859 


0.573 


-10.21 


-13.5 


-113.6 


63.0 


0.71 


0.80 


3.55 


4.03 


X Per 


51497.830 


0.577 


-10.22 


-12.9 


-115.1 


61.7 


0.68 


0.79 


3.52 


4.09 


X Per 


51817.899 


0.855 


-18.63 


17.5 


-42.7 


118.6 


2.29 


2.05 


3.17 


2.83 


X Per 


51818.900 


0.859 


-18.69 


16.0 


-43.1 


118.0 


2.30 


2.10 


3.14 


2.87 


X Per 


51819.887 


0.863 


-18.54 


15.7 


-43.0 


115.8 


2.34 


2.09 


3.13 


2.79 


X Per 


51820.918 


0.867 


-18.59 


15.7 


-43.5 


116.3 


2.31 


2.13 


3.13 


2.88 
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Star 


Date 


Orbital 


W x 


V r (W) 


V r (V) 


V r (R) 


W X (V) 1 




FWHM(V) 


FWHM(ii) 


Name 


(HJD-2,400,000) 


Phase 


(A) 


(km 


(km s- 1 ) 


(km s- 1 ) 


W X (R) 


V/R 


(A) 


(A) 
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3.06 


3.01 


A ir^cr 
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-18.97 


8.4 


-55.0 


105.2 


1.93 


1.85 


3.10 


2.97 


X Per 


51899.804 


0.183 


-18.93 


8.7 


-55.3 


104.1 


1.91 


1.84 


3.07 


2.96 


X Per 


51900.797 


0.187 


-18.83 


8.9 


-54.0 


105.5 


1.95 


1.89 


3.09 


2.99 


X Per 


51901.783 


0.191 


-19.26 


9.0 


-52.4 


107.2 


2.06 


1.97 


3.07 


2.93 


X Per 


53290.922 


0.740 


-17.71 


-9.7 


-80.7 


53.6 


0.66 


0.87 


2.65 


3.46 


X Per 


53290.926 


0.740 


-17.61 


-9.4 


-80.7 


56.8 


0.81 


0.92 


2.89 


3.28 


X Per 


53292.922 


0.748 


-18.12 


-9.3 


-78.3 


56.2 


0.78 


0.93 


2.74 


3.26 


X Per 


53292.928 


0.748 


-17.93 


-9.1 


-82.4 


51.7 


0.59 


0.83 


2.57 


3.64 


X Per 


54019.950 


0.653 


-23.81 


-6.4 


-8.9 








7.92 




X Per 


54020.939 


0.657 


-23.61 


-4.4 


-6.9 








7.78 




X Per 


54021.930 


0.661 


-24.38 


-4.5 


-7.2 








7.78 





Table 2 

Adopted Stellar and Disk Parameters 



Parameter 


HDE 245770 


X Per 


Rs (Re) ■■ 


15. 


6.5 


M s (M ) . 


20. 


15.5 


T s (K) .... 


28000 


29500 


i (deg) 


28.5 


31.5 


Pouter (Rs) 


15.9 


42.9 



Table 3 

IUE Radial Velocity Measurements 



Date 


Orbital 


V r 


(HJD-2,400,000) 


Phase 


(km s _1 ) 


43602.749 


0.034 


-5.3 


43705.644 


0.445 


3.3 


43712.142 


0.471 


20.4 


43784.061 


0.758 


6.2 


43796 362 


0.808 


7.5 


43850.204 


0.023 


4.3 


43885.683 


0.164 


5.8 


43885.749 


0.165 


-3.5 


43885.806 


0.165 


-5.6 


43947.719 


0.412 


1.6 


44156.035 


0.245 


-0.7 


44231.213 


0.545 


3.8 


44319.029 


0.896 


-3.2 


44319.087 


0.896 


-5.2 


44319.138 


0.896 


-0.3 


44319.187 


0.896 


-4.3 


44319.238 


0.897 


-4.3 


44319.279 


0.897 


-9.4 


44323.732 


0.915 


2.3 


44327.924 


0.931 


0.7 


44582.952 


0.950 


-0.6 


44582.991 


0.950 


4.2 


44888.132 


0.169 


-3.9 


48118.467 


0.075 


-4.4 


49593.258 


0.967 


5.8 


49593.325 


0.968 


2.1 


49621.183 


0.079 


-0.2 


49621.275 


0.079 


1.1 


49643.343 


0.168 


-4.0 


49647.365 


0.184 


-1.0 


49715.074 


0.454 


2.5 


49715.112 


0.454 


10.0 


49715.159 


0.454 


6.3 


49715.198 


0.455 


0.2 


49748.731 


0.589 


0.2 


49748.796 


0.589 


-0.9 


49775.778 


0.697 


8.2 


49775.824 


0.697 


7.0 


49956.299 


0.418 


2.9 


50000.094 


0.593 


-4.1 


50104.187 


0.009 


-15.1 


50124.212 


0.089 


1.2 


50144.055 


0.168 


1.0 
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Fig. 1. — A sequence of red spectra of HDE 245770 
from 1998 when the the disk had almost disap- 
peared. Each spectrum has its continuum aligned 
with the heliocentric Julian date of observation, 
and each is scaled in flux so that 10% of the con- 
tinuum equals 1 d of time. The bottom synthetic 
spectrum represents a model with T e fj = 28000 K, 
\ogg = 3.3, and Vsini = 230 km s _1 . The fea- 
tures present include the stellar Ha A6563 and He I 
A6678 lines and the interstellar features at 6613 
A and 6660 A. 
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Fig. 2.— The average spectra of HDE 245770 from 
each of the runs plotted with their continua off- 
set for clarity. Spectral intensity is scaled so that 
100% of the continuum equals four units of nor- 
malized flux. 
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Fig. 3. — The mean spectra of X Per from each 
run (separated in rectified intensity for clarity). 
Spectral intensity is scaled so that 100% of the 
continuum equals one unit of normalized flux. 
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Fig. 4. — The model relation between rescaled Ha 
emission equivalent width and the ratio of disk ra- 
dius (defined as the half maximum intensity radius 
of the integrated Ha flux) to the stellar radius. 



20 



2006 

+ 



0.4 r 

: o.2 1 



% + 4 



1000 2000 3000 4000 

DATE (HJD - 2,450,000) 



Fig. 5. — The time evolution of the X-ray flux 
(top), V (or J proxy) magnitude (middle), and 
disk radius (bottom) of HDE 245770. The disk 
radii are derived from the Ha equivalent widths 
and estimated V magnitudes (plus signs for data 
from Table 1 and crosses for values from earlier 
work). The dotted lines indicate disk radii associ- 
ated with the historic mean le vels o f emission iden- 
tified by iLvutv fc Zaitseval (|2000l ). The dashed 
line indicates the disk truncation radius for n = 5, 
and the dot-dashed line at the top marks a radius 
equal to the periastron separation. 
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Fig. 6.— X-ray light curves of HDE 245770 from 
RXTE/ ASM observations during quiescent dates 
formed by bi n ning w ith the orbital ephemeris from 
iFinger et alj ( 1994 ). The high energy band (top 
panel) appears to attain a maximum at a phase 
0.3 past periastron. 
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Fig. 7. — The time variations in X-ray flux (top 
panel: RXTE/ ASM counts binned in 25 d incre- 
ments), V magnitude (middle panel: AAVSO ob- 
servations binned in 25 d increments), and the ra- 
tio of disk to stellar radius (lower panel: based 
upon the rescaled Ha equivalent width). The plus 
symbols in the lower plot indicate radii derived 
from our equivalent width measurements, while 
the crosses are other published measurements. 



Fig. 8. — X-ray light curves of X Per from 
RXTE/ ASM observations during active dates 
form ed by binning with the o rbital cphcmeris 
from iDelgado-Marti et all (|200ll ). The high en- 
ergy band (top panel) appears to attain a maxi- 
mum at a phase w 0.25 past periastron. The verti- 
cal segments indicate the ±lcr standard deviation 
of the mean within each bin. 
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Fig. 9. — Radial velocities of X Per from IUE spec- 
troscopy plotted against the orbital ephemeris of 
the X-ray pulsar. The solid line shows the nominal 
best fit of the Be star systemic velocity and semi- 
amplitude (with all other parameters set by the 
pulsar orbit) while the dotted line shows a fit with 
the semiamplitude fixed at the 2a upper limit. 
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Fig. 10. — A mass plane diagram for X Per show- 
ing constraints on the Be star mass Mi and the 
neutron star mass Mi. The solid line shows the 
upper limit on mass ratio determined from the 
IUE radial velocities, and the dotted lines show 
the mass relations for three values of orbital incli- 
nation. The shaded region shows the most proba- 
ble range in masses for both components. 
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